Effects of gravitational confinement on bosonic asymmetric dark matter in stars 
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Considering the existence of old neutron stars puts strong limits on the dark matter/nucleon 
cross section for bosonic asymmetric dark matter. Key to these bounds is formation of a Bose- 
Einstein condensate (BEC) of the asymmetric dark matter particles. We consider the effects of 
the host neutron star's gravitational field on the BEC transition. We find this substantially shifts 
the transition temperature and so strengthens the bounds on cross section. In particular, for the 
well-motivated mass range of 5 — 15 GeV, we improve previous bounds by an order of magnitude. 



INTRODUCTION 

While the existence of dark matter is at present well es- 
tablished, the nature of this dark matter is still a matter 
of considerable speculation 1] . Constraints from preci- 
sion cosmology [2| and observations of the Bullet Clus- 
ter [3J] give a picture of the generic phenomenology that 
dark matter models must satisfy. Direct detection exper- 
iments have also placed constraints on the interactions of 
the dark matter with baryons. Recently, attempts have 
been made to use astrophysical observations to further 
constrain the nature of the dark matter. However, such 
constraints tend to be strongly model dependent. 

One particularly well-motivated class of dark matter 
models are known as asymmetric dark matter (ADM) 
models [4|, [5| . These models attempt to explain the sim- 
ilarity of dark matter and baryonic densities Udm ~ 
5f2baryon by postulating a mechanism for transferring 
the baryonic asymmetry of the standard model sector 
to the dark matter sector and visa versa. Such models 
favor a mass of the dark matter particles in the range of 
5—15 GeV [6]. The asymmetry in the dark sector implies 
that dark matter particles are stable and lack antiparti- 
cles with which to annihilate. 

This stability allows for interesting bounds to be de- 
rived from the existence of old neutron stars, particularly 
for bosonic dark matter @, [I]. The stability of ADM 
allows a neutron star to slowly accumulate a cloud of 
ADM particles by scattering them in such a way as to 
reduce their speed below the neutron star's escape speed. 
Once captured in the neutron star's gravitational field, 
repeated scattering between the star and the ADM par- 
ticles allows the ADM to thermalize with the star. Once 
a sufficient amount of dark matter has accumulated, the 
ADM cloud can become self-gravitating and collapse into 
a black hole. This black hole then consumes the host neu- 
tron star. Bosonic ADM can undergo a phase transition 
to a Bose-Einstein condensate (BEC) state. The BEC, 
having higher density than a thermal gas cloud, consid- 
erably increases the strength of the bounds obtained by 
this method. 

However, all analyses of this situation, to date, have 
treated the thermal ADM cloud as homogeneous. In this 



note, we consider the shortcomings of this approximation 
for a cloud of non-interacting bosons with a radius much 
smaller than that of the neutron star. Taking account 
of the effects of the confining gravitational potential im- 
proves the strength of the bounds for all masses up to 
15 GeV. The improvement is particularly strong, a fac- 
tor of an order of magnitude, in the well-motivated mass 
range of w 5 — 15 GeV. 



EFFECTS OF CONFINEMENT 

The well-known newtonian theory of gravitation in 
a uniform density sphere suggests that the ADM 
cloud within the neutron star should experience a har- 
monic potential energy given by V{r) — {2irGpm/2>)r 2 . 
Calculations^], subsequently verified by experiments 
with dilute alkali gases [10J, show the density of states 
for a gas trapped in a harmonic potential is markedly 
different from that of a homogeneous gas. Specifically, 
for a gas confined in an isotropic harmonic potential the 
density of states, g(e), is given by g(e) = e 2 /(2(huj) 3 ), 
where uj is the angular frequency of the harmonic oscilla- 
tor potential. This may be compared to the homogeneous 
gas result, g(e) = ^/e(2Tr(2m/ h 2 ) 3 ^ 2 V , where m is the 
dark matter particle mass and V is the volume occupied 
by the gas. 

For a fixed temperature, T, the critical number, N c , 
in the harmonic trap for which a BEC begins to form is 
N c = ((3)(kT/haj) 3 , with £(x) the Riemann zeta function 
and C(3) ~ 1.2. Using the newtonian potential above we 
find a critical number 



V 



sho 



C(3) 



kT 



hy/inGp/S 



Comparing this to the value, Vh om , found using the ho- 
mogeneous gas in a volume determined by gravitational 
binding in[7, 8] we see the same temperature scaling but 
a larger pref actor: 
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Using the newtonian potential for a relativistic sys- 
tem such as a neutron star, particularly near the core, 
misses important effects. We therefore derive the metric 
near the center of the star and from it derive an effective 
gravitational potential energy. To do so, we model the 
star as spherically symmetric with a metric of the form 
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where m(r) is the mass contained within a sphere of 
radial coordinate r. For the temperatures of interest, 
T « 10 5 — 10 7 K at the neutron star core, the thermal 
ADM cloud should be localized within a sphere of ra- 
dius roughly 1 m around the center of the neutron star. 
Thus it is reasonable to treat the density as uniform in 
the region of interest. This gives us m(r) — Airpr 3 /3, 
making g rr = 1 — 8nr 2 p/3. Using p = 1.5 X 10 15 g/cm 3 
shows that for r < 1 m, g rr — 1 < 10~ 8 . Thus, we 
can neglect the curvature of the spatial components of 
the metric and treat the radial coordinate as the proper 
length, simplifying the analysis. This allows us to treat 
the system via an effective newtoni an p otential, <fr. Using 
the Oppenheimer-Volkoff equation ll[ we find 
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where P is the pressure at the star's core. 

From this expression, we see that the newtonian analy- 
sis can be carried over by simply replacing p with p + 3P. 
The pressure depends sensitively on the neutron star 
equation of state. Recently, observations of x-ray emis- 
sion from neutron stars have progressed to the point 
where goo d a p proximations for this equation of state can 
be had[12l.ll3|. The core density used above corresponds 
to p « 0.9 GeV/fm 3 with a pressure of P w 0.3 GcV/fm 3 . 
The relativistic analysis therefore gives a gravitational 
potential roughly twice as strong as the newtonian anal- 
ysis. This further decreases the critical number for the 
condensation transition. Our final result for the critical 
number is thus 
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For the mass range 5—15 GeV the constraints on 
dark matter/nucleon scattering are dominated by the 
critical number. As such, using the corrected expression 



above strengthens the bounds in this range by roughly 
a factor of 8 as compared to those from ref [8j. Fur- 
ther refinements to the neutron star equation of state 
are anticipated in the near future, based on discrepancies 
between complimentary techniques for deriving it from 
observation [13l [lj] as well as the accumulations of fur- 
ther data. Bounds on dark matter cross-sections may 
be easily updated in the light of such new information, 
based on the above results. 
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